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Abstract 

The pseudo-critical temperature of the confinement-deconfinement transition and the phase transition surface are investigated by 
using the complex chemical potential. We can interpret the imaginary chemical potential as the Aharonov-Bohm phase, then the 
analogy of the topological order suggests that the Roberge-Weiss endpoint would define the pseudo-critical temperature. The 
I behavior of the Roberge-Weiss endpoint at small real quark chemical potential is investigated with the perturbative expansion. The 
expected QCD phase diagram at complex chemical potential is presented. 
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1. Introduction 

,__i Understanding the phase structure of Quantum Chromody- 
,1^ namics (QCD) is one of the most important and interesting sub- 
^-ijects in nuclear and elementary particle physics. The lattice 
^^QCD simulation is a powerful and gauge invariant method, but 
it has the sign problem at finite real chemical potential (//r), 

, and we cannot obtain reliable results at large /tr. Some meth¬ 
ods are proposed to circumvent the sign problem, see Ref. Ill] 
(N as an example, but those methods are limited in the i-Ir/T < 1 
region where T is temperature. Because of the sign problem, 
low energy effective models of QCD are extensively used to ex- 

_plore the QCD phase diagram. Effective models, however, have 

'sD strong ambiguities and thus quantitative predictions are impos- 
sible at present. Towards unification of lattice QCD simulations 
and effective model approaches, a new method so-called imag- 
inary chemical potential matching approach Hi is proposed 
lO recently. In this method, we use lattice QCD data obtained at 
finite imaginary chemical potential (pi) to constrain effective 
models. It is well known that the sign problem does not exist 
in the finite pi region and the region has information on the pR 
region; constrained models are reliable not only at finite pi but 
a also at finite pR. 

In addition to the imaginary chemical potential matching ap¬ 
proach, the concept of the imaginary chemical potential may 
be important to define the confinement-deconfinement transi¬ 
tion. In the imaginary time formalism where its time direction 
is compactified, the imaginary chemical potential can be inter¬ 
preted as the Aharonov-Bohm phase induced by U{1) flux in¬ 
sertions to the closed time loop. From this interpretation, we 
can determine the pseudo-critical temperature of the deconfine¬ 
ment transition from the Roberge-Weiss (RW) endpoint 0] with 
an analogy of the topological order fill as explained later. 

It is also interesting to investigate the phase structure at finite 
complex chemical potential, p - pR+ipi, where Pr + 0 and pi + 
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0. On the (T,pr) plane (pi = 0), the first order phase transition 
may exist, then we would have the critical point (CP) as the 
endpoint of the first order phase transition boundary. On the 
(T,pi) plane (pR = 0), the first order RW transition line exists, 
and the RW endpoint can either be the first or the second order. 
Some lattice QCD simulations fi suggest that the order 
of the RW endpoint is the first order at the physical pion mass. 
Then the RW endpoint is a triple point, with two other first order 
lines departing from the RW transition line. The first order lines 
appear around the heavy quark limit as well as around the chiral 
limit fii. These lines with small quark mass may have the 
chiral transition nature and are referred to as the chiral critical 
lines. One of the chiral critical lines extending in the pi —> 
0 direction should have an endpoint (chiral critical endpoint; 
CCE) before reaching the = 0 as long as the p — Q transition 
is crossover. Figure[T]summarizes our current expectation of the 
QCD phase diagram. Now we can raise a question. How does 
CCE behave at complex chemical potential ? Specifically, is 
CCE on the {T,pi) plane connected with CP on the {T,pr) plane 
at complex chemical potential, i.e. in the {T, pR, pi) space ? The 
topology of the phase diagram at complex chemical potential 
would tell us the relation between the deconfinement and the 
chiral transition. 

In this letter, first we briefly summarize properties of QCD 
at finite imaginary chemical potential, and propose a new def¬ 
inition of the pseudo-critical temperature of the deconfinement 
transition by the RW endpoint temperature {Trw)- Next, we in¬ 
vestigate the behavior of the RW endpoint at small pR by using 
the perturbative expansion. Finally, we present two scenarios of 
the QCD phase diagram at complex chemical potential based on 
the behavior of the RW endpoint at small pR and a symmetry 
argument. 

2. QCD with imaginary chemical potential 

At finite pi, QCD has a special periodicity so-called the RW 
periodicity 0. The RW periodicity can cause the first order 
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Figure 1: Schematic figure of our current expectation of the QCD phase di¬ 
agram at finite /jr and fi\, respectively. Solid lines represent the first order 
transition line. 


transition lines (RW transition lines) and their endpoints (RW 
endpoints). Those are predicted by using the strong coupling 
QCD and the perturbative one-loop effective potential with a 
background gauge field. The RW periodicity can be seen from 
the relation of the grand-canonical partition function Z |0]: 

Z(0)=Z(0+—), (1) 

where a dimensionless quark imaginary chemical potential is 
defined as 0 = fj.il T e R and k is any integer. If quarks are con¬ 
fined, physical states are classified by hadron degrees of free¬ 
dom unlike the deconfined phase. In the confined and decon- 
fined phases, origins of the RW periodicity are different: 

Confined phase — The origin is the dimensionless baryon 
chemical potential 30 in the form of exp(+3/0) in the parti¬ 
tion function. It can be easily seen, for example, from the 
strong coupling limit of the lattice QCD with the mean- 
field approximation MM, the chiral perturbation theory 
with the relativistic Virial expansion 11211 and that with the 
finite energy sum rule Q. 


Deconfined phase — The origin is the quark chemical poten¬ 
tial and the gauge field in the form of exp[+!(gA 4 / 7 ’ + 0 )]. 
It generates color non-singlet contributions in the parti¬ 
tion function. It can be understood from the perturbative 
one-loop effective potential with the background gauge 
field 114, and the RW periodicity is induced by Z 3 
images Then, we can find the nontrivial degeneracy of 
the free energy minima. 


The first order RW transition lines at 0 = 7T(2k + l)/3 are in¬ 
duced by Z 3 images. Its endpoint which is nothing but the RW 


endpoint should exist. The order of the RW endpoint is still un¬ 
der debate, but some lattice QCD simulations El [a Ist] suggest 
that the order seems to be first order around the physical pion 
mass. 


3. Imaginary chemical potential and Aharonov-Bohm 

phase 

In the imaginary time formalism where the temporal direc¬ 
tion is compactified, the t/(l) flux can be inserted to the closed 
loop of the temporal coordinate. Then, the imaginary chemical 
potential can be regarded as the Aharonov-Bohm phase 01611 . 
The Aharonov-Bohm effect has been discussed in the spatial 
loop with the flux insertion. In the imaginary time formalism, 
the temporal coordinate shares almost the same features as the 
spatial coordinates and thus we can use this interpretation. In 
this interpretation, we may use the discussion of the topolog¬ 
ical order j^]. Actual applications to zero temperature QCD 
was discussed in Ref. |@]. It should be noted that the Polyakov- 
loop (€>) is usually used to discuss the deconfinement transition. 
The Polyakov-loop can be obtained from a holonomy which is 
the gauge invariant integral along the temporal coordinate loop. 
The Polyakov-loop is the order parameter of the center symme¬ 
try breaking and it is related with the deconfinement transition 
in the infinite quark mass limit. The nontrivial degeneracy of 
the free energy minima in the deconfined phase may have some 
relation to the complex phase of the Polyakov-loop. It can be 
probed by the insertion of the flux as discussed below. 

In Ref. 101, the authors consider the torus at zero temper¬ 
ature, and introduce three adiabatic operations: (a) Insert the 
t/(l) flux to spatial closed loops, (b) exchange /-th and i + 1 -th 
quarks and (c) move a quark along loops. Commutation rela¬ 
tions of the operation (b) and (c) are described by the Braid 
group, and the Aharonov-Bohm effect determines the commu¬ 
tation relations of those with (a). If quarks are deconfined, op¬ 
erations become non-commutable because of the quark’s frac¬ 
tional charge. It is commutable if quarks are confined because 
physical states are described by hadron degrees of freedom with 
integer charges. Therefore, if there is only one vacuum in the 
deconfined phase, it is inconsistent with the non-commutability 
of the operations and thus vacuum degeneracy should exist. 

At finite temperature, the topological order cannot be well 
defined, because thermal states are constructed by a mixture of 
pure states with the Boltzmann factor, and we cannot operate 
(a), (b) and (c), adiabatically. Nevertheless, the RW periodicity 
shows significantly different behaviors in the confined and de¬ 
confined phases as already mentioned, and it is induced by the 
nontrivial appearance of the RW periodicity in the deconfined 
phase. This fact suggests that we can distinguish the confined 
and deconfined phases at /ii =0 from the non-trivial degen¬ 
eracy of the effective potential at 0 = n/3. We here consider 
the t/(l) flux insertion of 27r/3 to the temporal loop at zero fii. 
This corresponds to changing the t/(l) flux from -tt/3 to 7 r /3 at 
0 = n/3. These two states are free energy minima at 0 = 7 r /3 de¬ 
generated at r > Trw, while they belong to the same minimum 
at r < Trw. Thus the gluon configurations in states at 0 = n/3 
are essentially the same as those at 0 = 0. This degeneracy 
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seems similar to the vacuum degeneracy in zero T systems and 
the analogy can be found; the response of hidden local minima 
by the flux insertion via yUi tells us the non-trivial degeneracy of 
the free energy minima. Therefore, we propose that Trw is the 
pseudo-critical temperature of the deconflnement transition. 

Let us examine if Trw provides a deconflnement tempera¬ 
ture reasonably well. It should be noted that the present def¬ 
inition and the standard definition determined by using the 
Polyakov-loop are consistent in the infinite quark mass limit 
where the Polyakov-loop is the exact order-parameter of the 
deconflnement transition. Therefore, we can And the relation 
I'd = To = Trw where is the critical temperature deter¬ 
mined by the susceptibility of the Polyakov-loop and I'd is the 
deconflnement critical temperature. When the dynamical quark 
is taken into account, T® becomes the pseudo-critical tempera¬ 
ture. The upper bound of the pseudo-critical temperature can be 
determined by using the appearance of local minima of the ef¬ 
fective potential as found in the perturbative one-loop effective 
potential in the Re O - Im 4> plane. We call it Tz,. From the lat¬ 
tice QCD and effective model predictions, we can have the rela¬ 
tion T® < Trw < Tzj. While the Polyakov-loop is no longer the 
exact order parameter and thus the determination of I'd is not 
unique, Irw is uniquely determined. Thus, Irw which agrees 
with Id in the infinite quark mass limit is unambiguously de¬ 
termined in the lattice QCD and effective models and provides 
a reasonable value as Id with the dynamical quark. If we adopt 
Irw as Id, we lead to an implication that the deconflnement 
transition is the topological phase transition. 


4. Roberge-Weiss endpoint at complex chemical potential 

We now discuss the /iR-dependence of the pseudo-critical 
temperature of the deconflnement defined by Irw- We here give 
an model independent argument based on the perturbative ex¬ 
pansion of the effective potential in //r as a first step to investi¬ 
gate the /iR-dependence of the RW endpoint. It should be noted 
that non-perturbative model approaches have several difficul¬ 
ties. One of the promising effective models is the Polyakov- 
loop extended Nambu-Jona-Lasinio (PNJL) model 01711 . The 
PNJL model has the model sign problem at finite yUR 1 1811 . There 



directly used at finite complex chemical potential because we 
cannot maintain the RW periodicity and some other desirable 
properties of QCD. 

The effective potential at small yUR are expanded to yu^ order 
as 


T(I,yUR,yUl) 


- ^(1,0,yui) - (y)(In,(/tR,//i)U=o) 


1/PrC ^[I»?0^r,Pi)] | 
2 ^ 7 ’'^ dfip^lT i/iR=o 


0{(MR/Tf) 


( 2 ) 


where 


dna 


= I- 


dn„ 


dfj.RlTW=o d{ifiilT)W=o 


= -iT- 


driq I 

d6 l/iR=o' 


( 3 ) 


Equation Q is real, and the second term in r.h.s. of Eq. (|2]i 
should be pure imaginary. 

We here neglect the imaginary part of the effective potential. 
This assumption corresponds to the phase quenched approxi¬ 
mation. If the RW endpoint is the weak first order, the riq gap 
at Irw is small and the sign problem is mild, then the phase 
quenched approximation is justified. If the RW endpoint is the 
strong first order, the partition function at given (I,//) is dom¬ 
inated by one classical vacuum, then the role of the imaginary 
part, or the phase of the state, is minor. The phase quenched 
approximation cannot be applied to the second order RW end¬ 
point. Eortunately, the RW endpoint with the realistic quark 
mass is predicted to be the first order by lattice QCD simula¬ 
tions ItI [a Ist] and it would be possible to ignore the imaginary 
part of the effective potential. Effects of the imaginary part of 
the effective potential will be discussed elsewhere. 

Next, we consider the confinement and deconflnement con¬ 
figurations. We call the configuration at (Irw - e) confinement 
configuration which is labeled as C ^ where e is a infinitesi¬ 
mal positive value. Also, we call the configuration at (Irw + e) 
deconfinement configuration which is labeled as C+f. 

By comparing Re 'V with C-e to that with C+e in the e ^ 0 
limit, we can distinguish whether Irw decreases or increases if 
we can estimate the third term of Eq. (I2]i. Eor example in the 
lattice QCD and effective model calculations 1271 yj], we can 
see that Eq. (O is negative below Irw and it becomes moder¬ 
ate above Irw- The correction term makes Re 'V with C-e 
higher than that with C+^ because the correction term is then 
positive with C+f- This means that the first-order Irw decreases 
with increasing yUR at least in the small yUR region. This behavior 
is consistent with the decreasing behavior of the pseudo-critical 
temperature of deconflnement transition defined by using usual 
determinations; for example, see Refs IBSS- 


5. QCD phase diagram at complex chemical potential 

By taking into account our perturbative result and the sym¬ 
metry argument, we can sketch expected QCD phase diagrams 
at finite complex chemical potential. Phase diagrams expected 
from our present discussions are summarized in Eig.|2] Because 
of the RW periodicity, the phase structure should be periodic 
along the 0 axis. 

The RW transition line on the (1, 0) plane (//r = 0) may 
be topologically connected with the first order phase transition 
boundary on the (1,/tR) plane {0 = 0). Two of the first order 
transition lines starting from the RW endpoint have chiral tran¬ 
sition nature and are referred to as the chiral transition lines 101- 
Then it is not unreasonable to expect that the endpoint of the 
chiral critical line on the (1, 0) plane is connected with the QCD 
critical point on the (1,/tR) plane. In this case, the first order 
phase boundary on the (I,//r) plane forms a chiral transition 
surface in the (I,/tR,0) space, and connects the (I,yUR) plane 
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Figure 2: Two possible schematic QCD phase diagrams at finite complex 
chemical potential. The top (bottom) panel represents the correlated (uncor¬ 
related) case between the chiral and RW transition surfaces. 


and the {T, 0) plane. The RW transition line extends in the finite 
yUR region and forms an RW transition surface in the 0) 

space. The RW endpoint may reach T = 0 as shown in the top 
panel of Fig. |2]or it may deviate from the chiral transition sur¬ 
face at some temperature. There is the possibility that Trw line 
becomes smaller than the chiral critical surface at moderate /tr 
and finally becomes larger than the chiral transition surface.. It 
is deeply related with a strength of the correlation between the 
chiral transition surface and the RW transition surface. 

Another possibility is that the first order transition lines on 
the {T, 9) plane are topologically separated from the first order 
phase boundary on the iT,ju^) plane, as shown in the bottom 
panel of Fig. |2 Trw first decreases at small /tr, but does not 
goes across the chiral transition surface. In this case, the de¬ 
confinement transition represented by the RW endpoint on the 
(T, 6) plane has less relevance to the first order phase bound¬ 
ary, which would be the chiral transition, on the (r,yUR) plane. 
Therefore, we can call this possibility imcorrelated case and the 
former possibility correlated case. 

The QCD phase diagram at finite complex chemical pote ntial 
is related with the following subjects. (I) In Ref. 31], the 
authors use experimental data to construct the canonical parti¬ 
tion function. Then, ^rw at yUR = 0 is used to clarify realized 
temperatures in experiments through the Lee-Yang zero analy¬ 
sis 113211331] . In the analysis, ^rw at finite yUR should be related 
with zeros inside the unit circle on the complex quark fugacity 
plane. If so, there is the possibility that we can strictly deter¬ 
mine realized temperatures in experiments if we can systemati¬ 
cally understand the behavior of zeros. (II) The analytic contin¬ 
uation in QCD from the imaginary to the real chemical potential 
is usually performed on the plane. In the continuation, we 
may miss some information such as an inhomogeneous conden¬ 
sate S. The analytic continuation on the complex chemical 
potential plane may restore the information missing. 

6. Summary 

We have proposed that the Roberge-Weiss endpoint provides 
a reasonable deconfinement temperature. The imaginary chem¬ 
ical potential can be interpreted as the Aharonov-Bohm phase 
induced by t/(l) flux insertions and then the analogy of the 
topological order can be used. In the deconfined phase, we can 
find the degeneracy of the free energy by inserting the t/(l) 
flux, and these states belong to different minima. On the other 
hand, in the confined phase, we cannot find such non-trivial 
structure of the free energy. This suggests that we can distin¬ 
guish the confined and deconfined phases ?A0 - ^ijT = 0 from 
the non-trivial degeneracy of the effective potential. Then, Trw 
becomes the pseudo-critical temperature at // = 0. 

In this study, we have used the analogy of the topological 
order at finite T. The topological transition does not have the 
usual order parameter, but a relation with an entanglement and 
topological entropies has been discussed in the condensed mat¬ 
ter physics; for example, see Ref. 1^ . In QCD, such a relation 
is not clear, but it is an interesting direction to confirm the topo¬ 
logical nature of the transition. The calculation of the entan¬ 
glement entropy in QCD is possible by using the lattice QCD 
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simulation, and pioneering works were done in Ref 11361137 


in the quenched case. The calculation with the dynamical quark 
will provide us information to confirm present discussions. 

Using the perturbative expansion in terms of /tr, we have 
investigated the behavior of Trw at small yUR and then the de¬ 
creasing behavior of Trw is obtained. Based on these results, 
we presented two scenarios of the QCD phase diagram at finite 
complex chemical potential. First scenario which is the cor¬ 
related case has the strong correlation between the chiral and 
deconfinement transitions. The RW endpoint at finite yUR fi¬ 
nally reaches the 7 = 0. Then, the critical point can become 
more complex than the usual expectation since two more first 
order transition lines can be connected at the critical point. The 
second scenario is the uncorrelated case. The RW endpoint is 
separated from the chiral transition surface. 

Since the complex chemical potential is related with the Lee- 
Yang zero analysis and the analytic continuation to the finite jUr 
region, understanding the QCD phase structure at finite com¬ 
plex chemical potential may have impact to the beam energy 
scan program in heavy ion collider experiments, investigation 
of neutron star structures and so on. Our results are based on 
perturbative calculations and thus the non-perturbative checks 
should be done. 
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